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Fatal Effects of Charges on Stability of Black Holes 
in Lovelock Theory 

Tomohiro Takahashi 
Department of Physics, Kyoto University, Kyoto, 606-8502, Japan 

We study the stability of static spherically symmetric charged black holes under tensor 
type perturbations in Lovelock theory which is a natural higher dimensional generalization 
of Einstein theory. We derive the master equation for tensor type perturbations and present 
the criteria for the stability. Examining these criteria numerically, we show that a black hole 
T— I , with small charge has the instability if its mass is as small as extreme mass. Combined with 

T-H ■ our previous result that neutral black hole has no dynamical instability in odd dimensions, 

• this result suggests that charges have fatal effects on black holes in Lovelock theory. 

<N : 

■ §1. Introduction 
^ \ 

^ ' One of the most exciting predictions of the braneworld with large extra-dimensions 

■ is the possibility of higher dimensional black hole creation at the LHC.^* Then, it 
\ is important to study higher dimensional black holes. Especially, stability analysis 

f"^' of higher dimensional black holes is important because unstable black hole solutions 

qH' are not attractor solutions; that is, black holes with instability must not be created 

^ \ at the LHC. 

bX)' The stability of higher dimensional black holes in Einstein theory has been in- 

tensively studied. The most famous one is stability analysis of higher dimensional 
CN ! Schwarzschild black hole by Kodama and Ishibashi.'^ They have shown that higher 

' dimensional Schwarzschild black holes are stable for all type perturbations. They 

QQ ■ have also extended their works to charged black holes in higher dimensions and de- 

\ rived master equations.'^ Furthermore, for charged black holes, numerical studies 

' have also been done.® For the stability analysis of rotating black holes in higher di- 

mensions, a group theoretical method is developed,™' which is used for 5-dimensional 
rotating black holes with equal angular momentsP and also used to study the sta- 
bility of squashed black holes.'^HS The stability of a special class of rotating black 
holes in more than 5-dimensions is also studied .'^^^ 

As mentioned above, many studies have been done based on Einstein theory. It 
is no less important than these to examine the stability of higher dimensional black 
d • holes in more general gravitational theories. This is because Einstein theory is not 

the most general gravitational theory which contains terms only up to the second 
order derivatives of metric in the equations of motion. The most general theory 
in this sense is Lovelock theory.fi^ Furthermore, because black hole productions 
occur at the fundamental scale, Einstein theory is not reliable any more. In fact, it 
is known Einstein theory is only a low energy limit of string theory.^^^ In string 
theory, there are higher curvature corrections in addition to the Einstein-Hilbert 
term.nS Thus, it is natural to extend gravitational theory into those with higher 
curvature corrections in higher dimensions. Lovelock theory belongs to such class 
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of theories!^''^ Then, it is worthy to extend the stabihty analysis to this general 
Lovelock theory. 

In Lovelock theory, it is known that there exist static spherical symmetric black 
hole solutionJIS'lI^ (and also topological black hole solutions are founcP^). For this 
static spherical symmetric black holes, stability analyses under all type perturbations 
have been performedf^"^ It is shown that there exists the scalar mode instability 
in odd dimensions, the tensor mode instability in even dimensions and no instability 
under vector type perturbation in all dimensions. In 2nd order Lovelock theory, 
stability analysis is extended to black string.^^S They have shown that scalar pertur- 
bations have an exponentially decaying behavior under s-mode approximations. 

However, because black holes may be produced by proton-proton collision at the 
LHC, it is natural to suppose these black holes have charges. In Lovelock-Maxwell 
theory, there exist static spherical symmetric charged black hole solution, namely 
charged Lovelock black hole solution.'^Sli Then, it is important to extend the above 
discussion to this charged solution. 

In this paper, we study the linear stability of charged Lovelock black hole so- 
lutions. Our purpose is examining the response of instability to the charge and 
considering the black holes creation in Lovelock theory at LHC. In this paper, for 
the first step, we only concentrate on tensor type perturbation . Furthermore, we 
mainly consider the second order and third order Lovelock theory for simplicity, 
although it is important to consider higher Lovelock tevuisW^ 

The organization of this paper is as follows. In section [21 we review Lovelock 
theory and confirm the existence of asymptotically flat black hole solutions. In 
section [3l we consider tensor perturbations and denote the condition for stability 
of charged Lovelock black holes. The method using in this section is the same as 
our previous paper .f^D In section U we check the criteria presented in section [3] 
numerically. We mainly examine the second order and third order Lovelock theory 
in this section. The final section [5] is devoted to the conclusion. 

§2. Charged Lovelock Black Hole Solutions 

In this section, we review Lovelock theory and introduce charged black hole 
solutions. 

In Ref. [T3|) . the most general symmetric, divergence free rank (1,1) tensor is 
constructed out of a metric and its first and second derivatives. The corresponding 
Lagrangian can be constructed from m-th order Lovelock terms 
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where Rxn^'^ is the Riemann tensor in D-dimensions and 5^^^'^"'^'^^'^ is the General- 

pir^l"'Pm'^m 
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ized totally antisymmetric Kronecker delta defined by 

/ siii c ••• C \ 

J/^IA^^-Mp _ ^ ^ '^-1 ^-2 ••• 

Then, Lovelock Lagrangian in D-dimensions is defined by 

k 

m=0 

where we defined the maximum order k = [{D — l)/2] and 6m are arbitrary constants. 
Here, [z] represents the maximum integer satisfying [z] < z. Then, fixing k, n = D — 2 
satisfies n = 2k — 1 or n = 2k. Hereafter, we set bo = —2A, bi = ai = 1 and 
bm = cim/fn {m > 2), for convenience. 

Then, the action for Lovelock-Maxwell system is 



S 



r k "1 

fV^g -2A+Y. {^^rn} - I ^\f,.F^\ (2-2) 

L m=l J 



where F^^y = A^-y — A^-^ means field strength of electromagnetic field and A^ is 
vector potential. Taking variation of this action with respect to g^j^, we can derive 
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where is Lovelock tensor defined as 

k 



m=l 
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and 
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is energy momentum tensor of Maxwell field. Furthermore, by varying ()2-2p by Af^, 
we can get parts of Maxwell equations 



F^'-u = 0, 



(2-6) 



and the rests of Maxwell equations can be derived from the identity dF = 0, which 
means 



F 



[pu;X] 



0. 



(2-7) 



As shown in Ref 15), there exist static spherical symmetric solution of these 
equations. Let us consider the following metric 



ds' 



-f{r)dt^ + l/f{r)dr^ + r^j^dx'dx^ , 
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where jij means metric of S"" and n = D—2. Using this, Lovelock tensor is calculated 
as follows: 

1 

otherwise = , (2'9) 



where is defined as 

and W[ili\ is 



/(r) = 1 — r^V'(r) 



(2-10) 



f = [if V'-] + V' - 

[ am = am |np=r^('^ - P)} 

For Maxwell field, we assume spherical symmetric electric field, that is, 

F^'' = E{r), otherwise = . (2-11 

Substituting this ansatz for ()2-6p shows that E{r) satisfies 



ar(r"S(r)) = ^ E{r) = Vn(n - l)Q/r" , (2-12) 

where ^n(n — V)Q is constant of integral which means charge of black hole. The 
normalization factor \Jn{n — 1) is only for convenience. Note that ()2-lip and ()2-12p 
satisfies the identity ()2-7p . Hence, this is exact solution for Maxwell equations. The 
energy-momentum tensor for this solution is 

T* = = = 

* ^ 2 2r2'^ 

• ^ ^ ■ ^ n(n-l)Q^ ■ 
* 2 * 2r2'^ ^ ' 

otherwise = . (2-13) 
The results ()2-9p and ()2-13p leads Lovelock equation (|2-4p as 

1 f^n+iVFH'i^"- '^''^''"^^^^ 



2r"-i ^ 2 r2" ' 

otherwise = . (2T4) 

Note that these equations are not independent. In fact, a derivative of the first 
equation with respect to r leads the second equation. Therefore, it is sufficient to 
consider only the first equation. Integrating this equation leads 

= ^.-%^ M{r) . (2-15) 
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In ()2-15p . fi is constant of integral which is related to ADM mass as 

where we used a unit WnG = 1. We can get this result from the asymptotic behav- 
io^^ and also gain by a bzackground- independent formalism. ' 

In this paper, we want to concentrate on asymptotically flat, i.e. yl = 0, solutions 
with a positive ADM mass /i > because such black holes could be created at the 
LHC. We also assume that Lovelock coefficients satisfy 

am>0 , (2-17) 

for simplicity. Furthermore, for numerical calculation, we nondimensionalize all vari- 
ables. Our choice ai = 1 means only scale of length is not fixed. There exist many 
candidates, but we use ^/a2 for nondimensionalization in this paper. For example, 
radial r can be nondimensionalized as f = r j yfai- From /(r) = 1 — r'^tp{r) = 
1 — f'^{a2'ip{r)), 'ip{r) should be nondimensionalized as ip = a2ip and this nondimen- 
sionalization leads that eq. (|2-15p is expressed as 



where 



a 



m=3 



(Xm ~ _ — (n-l)/2 i/^i _ — (n-l)/2|^| /„ .. „n 

/x = a2 \Q\=Oi2 \Q\ (2-19) 



are nondimensionalized Lovelock coefficients, nondimensionlized mass parameter and 
nondimensionalized charge parameter respectivrely. After this, we miss tilde, that 
is, r means nondimensionlized radius f for example. 

It is easily seen that (j2-18p has many branches. However, under our assumptions 
(|2-17p . an asymptotic flat solution must exist. For example, in A; = 2 case, ()2-18p 
leads two branches 



Mr) - I + yi + - 2QV^ , . 

' " I -1 - Vl + 2/i/r"+i-2QVr2- ' ^'^^ 

In these two branches, the upper branch is asymptotic flat branch because this 
behaves as /(r) = 1 — r^'0(r) ~ 1 — fj,/r^~^ in asymptotic region. 

Against k = 2 case, it is difficult to solve ()2-18p in general cases. However, we 
can recognize that there must exist an asymptotic flat branch by using graphically 
method which is expressed in Fig{TJ In this figure, we fix r and seek the corresponding 
roots ^p, that is, tp{r). In detail, we set the both sides of ()2-18p as y and draw two 
curves y = W[ip] and y = M(r)(= const.) in ip — y plane. The cross points of these 
curves correspond tl^^r). In order to check the existence of asymptotic flat branch, 
we must take care two points. The first point is behavior of W[il'] in -0 > 0. Using 
our assumptions (|2-17p . it is easily seen that M^[0] = and W[ip] is monotonically 
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asymptotically AdS branch asymptotically flat branch 

Fig. 1. In this figure, the curve means y — W[ip] and the dotted hue is y = M{r). The cross points 
of these mean the roots of (l2T8p . Note that radial r is fixed in this figure so M(r) is constant. 



M(r) 




Fig. 2. This figure shows the behavior of M{r). M(r) is always negative in r < ro = (Q^ / ^y^^" 
M{r) is monotonically increasing in r < rmax ~ (2n(5^/((n + , M = Mmax = 

/x^^^ ( 2^qT^ ) at r = Tmax aud M(r) is monotonically decreasing in r > rmax- In 

asymptotic region r ^ oo, this function approaches 

increasing in -0 > 0. The second point is behavior of M{r), especially in asymptotic 
region. This function behaves as Figj2j From this figure, it is easily seen that M(r) 
is always positive in asymptotic region and M{r) — t- as r — t- oo. Combining these 
two points and Fig{Tl we can know that there must exist a cross point oi y = W[ip] 
and y = M{r) in -0 > as long as we consider asymptotic region; that is, a root of 
(|2-18p which satisfy ^/^(r) > must exist in asymptotic region. And this root also 
satisfies -0 — )■ as r — )• oo because M(r) — )■ in this limit. Considering eq. ()2T8p . 
■0 ~ means ~ /i/r"+-'^ — Q^/r^" and /(r) ~ 1 — /x/r""-*^ + Q^/r^""^ which means 
this root expresses asymptotic flat. Note that this result also reminds us that the 
constant of integration ^ is proportional to ADM mass. After this, we only consider 
this asymptotic flat branch. 

From FiglU this asymptotic flat branch behaves — t- V'max — ?• as r moves 
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(a) 



(b) 



Fig. 3. (a); W[tp] has no extreme value in tp < 0. (b): W[ip] has extreme values in tp < 0. Whether 
W[ip] has extreme value or not depends on Lovelock coefficients Um and k. 



oo —7- rmax — ^ ^0) ip becomes negative in r < tq because M(r) is negative in 
this region. When r becomes stih smaher, this branch runs into singularity where 
Kretschmann invariant R^upaR^^^^ diverges. This variable can be calculated as 

Rf.uXpR^"^' = f" + 2n^ + 2n(n - 1)^^^/^ 

from metric ansatz ()2-8p . Hence, there are singularities at r = and also exist where 
derivatives of /(r) diverge. Especially, for the latter, the first derivative of /(r) is 

/' = - = -2r^ - r'g^^ (-(n + 1)^ + 2n^) , (2-21) 

so Kretschmann invariant diverges at 9^ [•(/'] = 0, that is, singularities also exist 
where y = becomes extreme value. This means that position of singularity is 

different between Figl3]-(a) and Figl3]-(b). If W[ip] has no extreme values in ^ < 
(FigEKa)), there is a singularity at r = 0: if W[ip] has extreme values at tp{rs) < 
0(Figl3]-(b)), there is a singularity at r = (0 < < ro). In either case, there is a 
singularity in r < rg. 

From the viewpoint of cosmic censorship conjecture, there must exist horizons 
outside of singularities. The branch we concentrate on is static and asymptotic 
flat, so event horizon coincides with killing horizon. Then horizon radius rn is 
characterized as 

= firn) = 1 - rl^irn) = 1 - r|,ViT , (2-22) 

which means ipn must be positive. This result and behavior of tp show that ru, if 
exists, is larger than tq; that is, there exists no naked singularity if /{rn) = has 
solutions. Hence, what we have to examine is conditions for existence of horizons. 
Horizons are characterized by ()2-22p and -0/^ also satisfies 
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mini T 

Fig. 4. In this figure, we draw two fines y = B[ip] and y = in ip ~ y plane. The cross point 

means ipmini- Under our assumption (|2-17p . y = B[ip] is monotonic function which diverges 
near t/j ~ and approaches as oo. 



Therefore horizons are determined by two equations 



H 



^ H 



' H 



(2-23) 



Substituting the first equation for the second equation leads 



(2-24) 



The solutions of this equation correspond horizon, so the behavior of y = 
determines the condition for existence of horizons. The derivative of Alip] with 
respect to ijj is 



n 



lg2^{n-3)/2 _ !i±l^-(n+3)/2^[^] ^ ^-(n+l)/2 « (2-25) 



2 " 2 
and solving 9^74[-0] = with respect to leads 



n — 1 
2 



n + 1 



n - 
= B\ 



1 



2 

n+1 



-n+l 



1 1 + 



n + 1 



ni=3 



1 ] ^ y- ^ 



n+l 



(2-26) 



Note that n = 2k — 1 ox n = 2k. Then powers of all terms in -B['i/'] are negative 
and their coefficients are positive. Therefore, y = B[ip] is monotonically decreasing 
function like FigHl From this figure, there must exist a root of (|2-26|) . We call this 
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i^mini because A[ip] becomes extreme minimum at ipmini which can be seen from 
()2-25p . Then, A[ip] behaves as FigJSl In this figure, we solve eq. ()2-24p graphically. 
The cross points of y = ^ and y = ^[V'] mean horizons. Prom this figure, it can be 
easily seen that there exist horizons if /U ^ /^ex where /Xe^ is defined as 

/^e. = Q^^LTJ^/' + ^d™^'^/ V[^_.] . (2.27) 

Note that ^l^out and ipm in Figl2] means outer horizon and inner horizon respectively 
because ipn corresponds horizon radius as xpn = ^/'''h- 

Finally, for example, we examine /iex in 5-dimensions (n = 3 and k = 2 case). 
From eq. (|2-26p . ipmini in 5-dimensions is 

2 1 1 

SO substituting this for ()2-27p leads 

/iex = 2|Q| + ^ . (2-28) 

§3. Condition for Stability under Tensor Type Perturbations 

In this section, we consider conditions for stability under tensor type perturba- 
tions. Note that we only consider /i > which means there exists horizon and we 
only consider the perturbation outside of the outer horizon. 

The background metric (j2-8p has n dimensional spherical symmetry, so pertur- 
bations are decomposed into scalar type, vector type and tensor type. In this paper, 
we concentrate on tensor type perturbations which are characterized as 

^gtt = Sgtr = Sgu = ^grr = ^gn = 

6g^j = r\{r)e'^'hl , (3-1) 
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where x means master variable and hjj is tensor harmonics. Tensor harmonics hfj 
satisfies f^hfj = 0, hf-\^ = and = -^thfj where -ft = i{i + n - 1) - 2 

{I = 2,3,4, •••). Here | means covariant derivative with respect to 5" metric ^ij. 
Note that there is no tensorial perturbation for Maxwell field A^, which means that 
the first order perturbation of energy momentum tensor 6T^^ = 0. Therefore, the 
first order perturbation of EOM is SQ^'^ = and this can be calculated as followsj^ 

n / // / 9 // 2fT' ' '\ > e(i + n-l)f „ 9 , , , 

-fTx - [fT +^ + //r ]x + ^j——^T x = UJ^T X . {3-2) 

\ r I [n — 2)r 

Here T(r) is defined as 

r(r) = r"-^d^W[ij] . (3-3) 

Now we will present the condition for the stability of the solutions we are con- 
sidering in this paper. 

As we will soon see, the master equation (j3-2|) can be transformed into the 
Schrodinger form. To do this, we have to impose the condition 

r'(r) > (for r>rH) . (3-4) 

In fact, this is necessary for the linear analysis to be applicable. In the case that 
there exists rg such that T (r^) = and rg > rn, we encounter a singularity. Using 
approximations T (r) ^ T {rg){r — rg) = T {rg)y, f{r) = f{rg) and r = rg, ()3-2p 
approximately becomes 

y^ + ^ + cx = o. (3-5) 

o-y ay 

This shows that near r = rg, x behaves as x ~ ci + C2logy , where ci and C2 
are constants of integration. Hence, the solution is singular at y = for generic 
perturbations. The similar situation occurs even in cosmology with higher derivative 
terms.l2S}'E3 In those cases, this kind of singularity alludes to ghosts. Indeed, if 
there is a region T (r) < outside the horizon, the kinetic term of perturbations has 
a wrong sign. Hereafter, we call this the ghost instability. 



When the condition (j3-4p is fulfilled, introducing a new variable !P'(r) = xi^)"!" (r) 
and switching to the coordinate r*, defined by dr* /dr = 1//, we can rewrite Eq. ()3-2p 
as 

-^2+^Mr*m = ^^^ , (3-6) 

where 

l{i + n-l)fd\nT' , 1 rdf.d j-\ 
^'^'^ = in - 2)r dT + ^^Tr Vd^'^^ ) ^^'^^ 



is an effective potential. 
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For discussing the stability, the "S-deformation" approach is usefulP''® Let 
us define the operator 

dr* 



(3-8) 



acting on smooth functions defined on / = (r|^, c«). Then, (j3-6p is the eigenequation 
and is eigenvalue of T-L. We also define the inner products as 

{^1,^2) = j <f*i'^2dr* . (3-9) 

In this case, for any smooth function ip with compact support in /, we can find 
a smooth function S such that 

i<f,Hip)= J^i\D^\^ + V\^\^)dr*, (3-10) 



where we have defined 



« = A + S, v- = y, + /f-5^ (3.11) 



Following Ref. [T8|) , we choose S to be 



S = -f-^ln{rVr) . (3-12) 



Then, we obtain the formula 



{ip,n^)= I \Dip{'dr*+eii + n-l) I -^^^^}^dr . (3-13) 



rH 



(n — 2)r dr 



Here, the point is that the second term in (|3-13p includes a factor i{i + n — 1) > 0, 
but T does not include i. Hence, by taking a sufficiently large £, we can always 
make the second term dominant. 

Now, let us show that the sign of d In T /dr determines the stability. If d In T /dr > 
on /, the solution ()2-8p is stable. This can be understood as follows. Note that 
i{i + n — 1) > 0, then we have V > for this case. That means {<p,T-Lip) > for ar- 
bitrary (/3 if dlnT /dr > on /. We choose, for example, ip as the lowest eigenstate, 
then we can conclude that the lowest eigenvalue Wg is positive. Thus, we proved the 
stability. The other way around, if dlnT /dr < at some point in I, the solution is 
unstable. To prove this, the inequality 

is useful. This inequality is correct for arbitrary ip, which is not necessarily an eigen- 
function of Ti, as far as the left hand side of Eq. (j3-14p make sense. If dlnT' /dr < 
at some point in /, we can find (p such that 

1^1' '""^-*<0. (3.15) 



(n — 2)r dr 
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In this case, ()3-13p is negative for sufficiently large i. Then, the inequality ()3-14p 
implies loq < and the solution has unstable modes. Thus, we can conclude that 
the solution is stable if and only if din T /dr > on I. 

From the above logic, if din T /dr has a negative region, negative cj^ states exist. 
Therefore, this instability is dynamical. Then, we call this as dynamical instability 
in order to distinguish this from the ghost instability which is caused by negativity 
of r'(r). 

We want to summarize this section. If T has negative region outside the outer 
horizon r > Vout, this solution has the ghost instability. Even if T is always posi- 
tive, this solution has the dynamical instability if T has a negative region outside 
the outer horizon. Therefore, charged Lovelock black holes are stable under tensor 
perturbations if and only if T and T are always positive outside the outer horizon. 

Note that T is calculated as T = r*^" ^1 + 2M(r) in A; = 2 case from ([2^ . 
This leads 



T' 



^l + 2M(r) 



(n-l) + (n-3)-^ + 2 



„2n 



>o. 



(3-16) 



which means there exist no ghost instabilities in 5 and 6 dimensions. 



§4. Numerical Results 



In this section, we examine "whether T' has negative region or not" and "whether 
T" has negative region or not" numerically. The detail is as follows. Note that 
Afi,dfj.,d\Q\ and Qend are parameters for numerical calculations. 

Firstly, we fix \Q\. Using this \Q\, extreme mass Hex can be determined from 
(|2-27p . Then, we change fi from ^Ugx to yUex + by d/i. For each fi, we can gain ■ip{r) 
with (IQI, /x) numerically from (j2T8p and determine T'{r) and T"{r). Then, checking 
whether these functions have negative region or not, we determine, for example, the 
border between stable and unstable which we call ^max- Then, we change \Q\ by 
d\Q\ and do the same calculation. We repeat this calculation until | Qend I- 

Note that we examine the region rout ^ r < rout + 10 and the mesh size dr is 
dr = 1.0 X 10-^ 

4.1. 5- dimensional Case 

Numerical results of 5-dimensional case are FigEJand Figd The former is Q — /i 
diagram near \Q\ ~ and the latter is that of \Q\ ~ 3. It can be seen that black holes 
with \Q\ = are stable, which agrees with our previous results.'^D When black holes 
are a little charged up, however, there exists an unstable region near extreme mass 
(Figi6|) and this region vanishes in \Q\ > 3 (FiglT]). As we have already discussed, 
there is no ghost instability in 5-dimensions. Hence, these figures show that nearly 
extreme black holes are unstable if < \Q\ < 3. 

4.2. 6-dimensional Case 

FigJS] and Figl9] are the numerical results in 6-dimensions. The former is Q — 
diagram near \Q\ = and the latter is that of \Q\ ~ 3.28. When black holes are 
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5D near |Q|=3 





1 




M-max 

extreme mass 






stable 










/^3,6.5) 




1 




no horizon 




2.8 


2.9 


3 3.1 3.2 



|Q| 



Fig. 7. Q — jj. diagram near \Q\ = 3 in 5-dimensions. This is calculated with — 0.05 and 
d/i — d\Q\ = 1.0 X 10~^. In this figure, the difi'erence between /imax and extreme mass is 
O(10~^) where \Q\ < 3 and there is no instability where \Q\ > 3. 



neutral, there is unstable region in < ^ < 0.27; this agrees with our previous 
results.!^ When black holes are a little charged up, there also exist an unstable 
region (FigjS]). However, this region vanishes in \Q\ > 3.282 (FigE]). As we have 
already discussed, there is no ghost instability in 6-dimensions. Therefore, it is 
shown numerically that nearly extreme black hole is unstable when its charge satisfies 
< IQI < 3.282. 
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6D near |Q|=0 



0.2709 
i 0.25 





eta hlo 






unstable 








no horizon 






extreme mass 






1 i^max 1 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 
|Q| 

Fig. 8. Q — /i diagram near \Q\ = in 6-dimensions. This is calculated with — 0.05 and 
dfi = 1.0 X lO""*. And d\Q\ is d\Q\ = 1.0 x 10"" when \Q\ < 10"^ and d\Q\ = 1.0 x 10"^ in 
101 > 10-^ 




Fig. 9. Q — fj, diagram near \Q\ — 3.28 in 6-dimensions. This is calculated with Ajj. = 0.01 and 
d/x = d\Q\ — 1.0 X 10""^. In this figure, the difference between fimax and extreme mass is O(10~"^) 
in |Q| < 3.282 and there is no instability in \Q\ > 3.282. 



4.3. 7- dimensional case 

The previous works tell that neutral black holes do not have ghost instability if 
C3 < 0.25 and have this instability if C3 > 0.25. Then, in this paper, we concentrate 
on C3 = 0.2 case and C3 = 0.3 case. 

4.3.1. C3 = 0.2 case 

FigJTO] and FigJTT] are the numerical results when C3 = 0.2 in 7-dimension. The 
former \s Q — ^ diagram near \Q\ = and the latter is that of \Q\ ~ 4.7. In 
FigllOt there are no unstable region and no ghost region, which are consistent with 
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our previous results.'^ However, when a little charged up, black hole has dynimical 
instability if its mass is as small as extreme mass. And this instability vanishes when 
IQI > 4.695 (Fig llip . Note that ghost region cannot be detected by our numerical 
calculation with < |Q| < 5, A^=l and dfi = d\Q\ = 0.01. Therefore, same as 
5-dimensional case, these result means that there exists an unstable region near 
extreme mass in < \Q\ < 4.695. 




Fig. 10. Q — n diagram near \Q\ = when C3 = 0.2 in 7-dimensions. This is calculated with 
Afi = 0.15 and dfi = 1.0 x lO""*. And d\Q\ is d\Q\ = 1.0 x 10"" when \Q\ < 10"^ and 
d\Q\ = 1.0 X 10"^ in \Q\ > 10"^ 



7D C3=0.2 near |Q|=4.695 



11.2 
11 
10.8 
10.6 
10.4 
10.2 
10 
9.8 



stable 


III. 




(4.695,10.5249) 




no horizon 




extreme mass 




1 1 ^ma)^ 



4.4 



4.5 



4.6 



4.7 

|Q| 



4.8 



4.9 



Fig. 11. Q — fi diagram near \Q\ — 4.7 when C3 — 0.2 in 7-dimensions. This is calculated with 
= 1.0 X 10"^, dfi = 1.0 X 10"* and d\Q\ = 1.0 x 10"^. In this figure, the difference between 
fJ-max and extreme mass is 0(10"*) in \Q\ < 4.695 and there is no instability in \Q\ > 4.695. 
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4.3.2. C3 = 0.3 case 

Fig{T2] and FigflBl are numerical results when C3 = 0.3 in 7-dimensions. The 
former is Q — fi diagram near \Q\ = and the latter is that of \Q\ ~ 5.4. In Fig ll2l 
there exists a ghost region if \Q\ = and 0.1 < /U < 0.15. This result agrees with the 
previous work.'^ As black hoes are charged up, this ghost region diminishes and 
unstable region appears. And this unstable region vanishes if \Q\ > 5.422 (Fig JTS]) . 

7D 03=0.3 near |Q|=0 

0.22 

0.2 

0.18 

0.16 
i 0.15 
0.14 

0.12 

0.1 



0.08 

0.002 0.004 0.006 0.008 0.01 

|Q| 

Fig. 12. Q — /i diagram near \Q\ — when C3 = 0.3 in 7-dimensions. fimax is calculated with Afi — 
0.1 and dfi — 1.0 x 10^*. Ghost and extreme mass are done with A/j, = 0.1, d^, = 1.0 x 10"'* 
and d\Q\ = 1.0 x 10"^ This fi gure shows that there exists a ghost region if \Q\ < 6.0 x 10 
and there exists an unstable region if IQ] > 4.9 x 10~*. 



7D 03=0.3 near |Q|=5.422 
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12.6 
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11.6 

5.2 5.3 5.4 5.5 5.6 5.7 5.8 

|Q| 

Fig. 13. Q — /i diagram near \Q\ — 5.4 when C3 = 0.3in 7-dimensions. This is calculated with 
Afi = 1.0 X 10"^, dfj, = 1.0 X 10"* and d\Q\ = 1.0 x 10"^. In this figure, the difference between 
fJ-max and extreme mass is 0(10^*) in \Q\ < 5.422 and there is no instability in \Q\ > 5.422. 
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4.4. 8- dimensional Case 

If black holes are neutral, in 8-dimensions, the previous works show that there 
is no ghost instability when C3 < 5.92 and there exists if C3 > 5.92. Therefore, we 
concentrate on C3 = 5.9 case and C3 = 6.0 case in this paper. 

4.4.1. C3 = 5.9 case 

FigHHand Fig(T5]are the numerical results when C3 = 5.9 in 8-dimension. The 
former is Q — fi diagram near \Q\ =0 and the latter is that of |(5| ~ 78.4. In 
FigHH black holes with Q\ = and < n < 4.513 are unstable, which agrees 
with our previous works.'^D When a little charged up, black hole has also instability 
if its mass is as small as extreme mass. However this instability vanishes when 
IQI > 78.386 (Fig JTS]) . Note that ghost cannot be found in our numerical calculation 
with < IQI < 80, Z\/U=5 and dfi = d\Q\ = 0.1. Therefore, same as 6-dimensional 
case, these result means that there exists an unstable region near extreme mass in 
< IQI < 78.386. 



8D C3=5.9 near |Q|=0 




2 - 

1 no horizon 



I 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

|Q| 

Fig. 14. Q — /i diagram near \Q\ = when C3 = 5.9 in S-dimensions. This is calculated with Ajj, = 5 
and dfi = d\Q\ = 1.0 x 10"^. 
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8D C3=5.9 near |Q|=78.386 



stable 



*^^(78.386, 168.244) 
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i 168.5 
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167.5 
167 

78 78.1 78.2 78.3 78.4 78.5 78.6 78.7 78.8 78.9 79 
|Q| 

Fig. 15. Q ~ n diagram near \Q\ — 78.4 when ca — 5.9 in 8-dimensions. This is calculated with 
Afi = 0.01 and dfi = d\Q\ = 1.0 x 10~^. In this fi gure, the difference between Umax and extreme 
mass is O(10~^) in \Q\ < 78.386 and there is no instability in IQj > 78.386. 



no horizon 



extreme mass 
^rnax 



4.4.2. C3 = 6 case 

FigHSland Figll7lare numerical results when C3 = 6 in 8-dimensions. The former 
is Q — fj, diagram near \Q\ = and the latter is that of IQI ~ 79.8. In FigIT6l there 
exists a ghost region if \Q\ = and < ;U < 4.888. This result agrees with our 
previous work.l^ As black hoes are charged up, this ghost region diminishes and 
unstable region appears. And this unstable region vanishes if \Q\ > 79.797(Fig JT7]) . 



8D 03=6.0 near |Q|=0 
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Fig. 16. Q — j-i diagram near |Q| = when C3 = 6 in 8-dimensions. fimax is calculated with Ajj, — 5.0 



d\Q\ = 1.0 X 10 . Ghost and extreme mass are done with = 5.0, dfi = 1.0 x 10 
= 1.0 X 10~*. This figure shows that there exists a ghost region if \Q\ < 5.2 x 10" 



and dfj, - 
and d\Q\ 

and there exists an unstable region if \Q\ > 8.8 x 10~ 
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8D C3=6.0 near |Q|=79.797 

173 p 
172.5 - 
172 - 

171.5 - ..- ■ no horizon 

(79.797,171.236) 

171 - 

extreme mass 

170.5 I 1 1 1 1 1 1 1 1 1 1 

79.5 79.6 79.7 79.8 79.9 80 80.1 80.2 80.3 80.4 80.5 

|Q| 

Fig. 17. Q ~ n diagram near \Q\ — 79.8 when C3 = 6in 8-dimensions. This is calculated with 
Afi = 0.01 and dfi = d\Q\ = 1.0 x 10"^. In this fi gure, the difference between Umax and extreme 
mass is O(10~^) in \Q\ < 79.797 and there is no instability in \Q\ > 79.797. 

4.5. Summarizing Numerical Results and Consideration 

In this subsection, we summarize features of these numerical results and consider 
black hoes at the LHC briefly. After this, we call both "ghost instability" and 
"dynamical instability" as only "instability" . 

In neutral case, there exists instability near // = if the spacetime is even 
dimensions and no instability in odd dimensions with small C3, which agree with 
the previous results. However, a little charged up, black hole suffers instability if its 
mass is near extreme mass in both even and odd dimensions, which is independent 
of Lovelock coefficient C3. Then, we can denote that charge is crucial for Lovelock 
black hole in odd dimensions. Another feature is existence of Qmax such that there 
is unstable region near extreme mass in \Q\ < Qmax and it vanishes in \Q\ > Qmax- 
This shows that extreme Lovelock black hole must not exist if its charge is less than 

Qmax ■ 

Using these results, we now consider charged Lovelock black holes at the LHC 
briefly. Because unstable black holes are never created, these events must occur on 
stable region in Q — /i diagram. So the numerical results show that mass has lower 
bound Umax in black hole creation \Q\ < Qmax- Or if we could observe a zero 
temperature black hole at the LHC, such black hole should has more charge than 

Qmax ■ 

After creation, black hole loses its charge quickly by Schwinger-type mecha- 
nisnpSl' and ran into evaporating phase. It is known that charge of black hole fluc- 
tuates in this phase.l^ As we have mentioned, black hole has instability if it has 
non-zero charge. Then, if dimension is even, black holes always run into unstable 
region; even if dimension is odd, black holes could become unstable by these fluc- 
tuations in this phase. That is, independent of dimension, black hole has always 
possibility to become unstable due to fluctuating charge in Lovelock theory while 



stable 
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they are evaporating. 

§5. Conclusion 

We have studied the stabihty of static charged black holes under tensor pertur- 
bations in the second order and the third order Lovelock theory. We have derived 
master equation and presented the criteria for stability. Examining these criteria 
numerically, we have shown that there exists an unstable region near extreme mass 
in both even and odd dimensions. In odd dimensions, especially, we have already 
shown that there is no dynamical instability when black holes are neutral. Therefore, 
charges have proved fatal to Lovelock black holes especially in odd dimensions. 

We have also shown that this instability is stronger as i becomes larger. And 
we have also pointed out that black holes have possibility to ran into the unstable 
region due to fluctuating charges in evaporating phase at the LHC. These suggest 
that black holes in Lovelock theory should be doomed to ruin at the LHC. 

One of future works is extension to vector and scalar type perturbations. For 
vector type perturbations, especially, we have already shown that black holes are 
stable in neutral case. Therefore, by checking vector type perturbations, it may 
become clearer that charges are perils for Lovelock black holes at the LHC. 

Extension to still higher dimensional case is also important. This is because 
string theory is consistently formulated only in 10-dimensions. Although the criterion 
for stability we have shown in this paper are true for any order Lovelock theory, this 
work is a little bother only because other parameters Cm are added. 

It is interesting to examine the nature of this instability. The shortest way may 
be search for the meaning of the function T(r) because this governs the dynamical 
stability of Lovelock black holes. Therefore, if T(r) has thermodynamical mean- 
ing, the relation between thermodynamicaP^'I^Sl' and dynamical instability might 
be revealed. Another way is considering physical properties. In Lovelock theory, 
gravitational collapses have also been studied.!^ They have shown that apparent 
horizon is less liable to appear in Lovelock theory, which means attractive force be- 
comes weaker; that is, higher curvature corrections are repulsive effectively. And 
this effective repulsive force should be dominant in large i mode. Then, origin of in- 
stability may be that repulsive force defeats attractive one, which is consistent with 
the results of this paper. This repulsive instability might exist in Einstein theory 
with galileons, too.'^S' 
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